NYU-TH-94/10/03 
October 1994 

The heavy quark decomposition of the S-matrix 
and its relation to the pinch technique. 

Joannis Papavassiliou, Kostas Pbilippides 
and 

Martin Schaden 

Department of Physics, New York University, 4 Washington Place, 
New York, NY 10003, USA. 

ABSTRACT 

We propose a decomposition of the S-matrix into individually gauge invariant sub-ampli- 
tudes, which are kinematically akin to propagators, vertices, boxes, etc. This decompsition 
is obtained by considering limits of the S-matrix when some or all of the external particles 
have masses larger than any other physical scale. We show at the one-loop level that 
the effective gluon self-energy so defined is physically equivalent to the corresponding 
gauge independent self-energy obtained in the framework of the pinch technique. The 
generalization of this procedure to arbitrary gluonic n-point functions is briefly discussed. 
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The pinch technique (PT) [1] is an algorithm that allows the construction of modified 
gauge invariant (g.i.) n-point functions, through the order by order rearrangement of the 
Feynman graphs contributing to a certain physical, and therefore ostensibly g.i. amplitude 
(such as an S-matrix element, or a Wilson loop). Even though the most recent applications 
of the PT are inspired by Standard Model phenomenology [2-7], it was originally intro- 
duced in the context of QCD, as a first step toward the construction of Schwinger-Dyson 
equations, which would respect the crucial property of gauge invariance, even in their one- 
loop dressed truncated version [8-9]. The simplest example that demonstrates how the 
PT works is the gluon two point function (propagator). Consider the .S-matrix element 
T for an elastic scattering process such as qiq 2 — > ?i?2, where 51,52 are two on-shell test 
quarks with masses m\ and 1712. To any order in perturbation theory T is independent 
of the gauge fixing parameter £. On the other hand, as an explicit calculation shows, the 
conventionally defined proper self-energy [collectively depicted in graph 1(a)] depends on 
£. At the one loop level this dependence is canceled by contributions from other graphs, 
such as 1(b), 1(c), 1(d), and 1(e) which, at first glance, do not seem to be propagator-like. 
That this cancellation must occur and can be employed to define a g.i. self-energy, is 
evident from the decomposition: 

T(s,t,m 1 ,m 2 ) = T (t,£) + T 1 (t,m 1 ,{ i ) + T 2 (t,m 2 ,C) + T 3 (s,t,m 1 ,m 2 ,C) , (1) 

where the function To(t,£) depends kinematically only on the Mandelstam variable t = 
— (pi — Pi) 2 = —q 2 , and not on s = (p\ + p 2 ) 2 or on the external masses. Typically, self- 
energy, vertex, and box diagrams contribute to To, T\ and T 2} and T3, respectively. Such 
contributions are £ dependent, in general. However, as the sum T(s,t,mi,m 2 ) is g.i., it is 
easy to show that Eq. (1) can be recast in the form 

T(s,t,mi,m 2 ) = f (t) + fi(t,mi) + f 2 (t,m 2 ) + f 3 (s,t,mi,m 2 ) , (2) 

where the Ti (i = 0, 1, 2, 3) are individually ^-independent. The propagator-like parts 1(f), 
1(g), 1(h), and l(i), stemming from graphs 1(b), 1(c), 1(d), and 1(e), respectively, enforce 
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the gauge independence of To(t), and are called "pinch parts". They emerge every time a 
gluon propagator or an elementary three-gluon vertex contributes a longitudinal to the 
original graph's numerator. The action of such a term is to trigger an elementary Ward 
identity of the form f = — m) — (J> — m) when it gets contracted with a 7 matrix. The 

first term removes the internal fermion propagator (e.g. it produces a "pinch"), whereas 
the second vanishes on shell. From the g.i. functions Ti (i = 1,2,3) one may now extract 
a g.i. effective gluon (G) self-energy fl^^g), g.i. Gqiqi vertices and a g.i. box B, in 

the following way: 

T = g 2 u 1 j^u 1 [(\)tl fliy (q)(\)]u 2 J 1 'u2 
q^ q^ 

r (3) 

T 2 = g 2 u 1 j^u 1 ( — )u 2 fi 2) U2 
q z 

f 3 =B 

where Ui are the external spinors, and g is the gauge coupling. Since all hatted quantities 
in the above formula are g.i., their explicit form may be calculated using any value of the 
gauge-fixing parameter £, as long as one properly identifies and allots all relevant pinch 
contributions. The choice £ = 1 simplifies the calculations significantly, since it eliminates 
the longitudinal part of the gluon propagator. Therefore, for £ = 1 the pinch contributions 
originate only from momenta carried by the elementary three-gluon vertex of graph 1(b) 
(and its mirror graph, not shown). The one-loop expressions of II ^^(q) and f ^} are given 
by [2], [9]: 

n F ( 9 ) = nf)( 9 )+vn p (9) (4) 

with t^y = q 2 g^ - q^q„, and 

nP «> = - 2l ^7„F<F 

t|-^-)[CW - — v 

l0TV z fi 



2c a(^)[CuV-H^-) + 2] 



(5) 
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where Cjjv = ~ — 7 + ln(47r), e = n — 4 ,7 = 0.577... is the Euler constant, j n = j (f^y™ 



is the dimensionally regularized loop integral. Similarly, for the vertex we have: 
[f,P =ig 2 

+ 1 



t c a , f l p S t (p t - k)j a T^ atl Ca f j p S z (p z - k)j tl S z ( Pz - k)j p 

(v) / ,.2//. , ^2 + (l7- c f) 



2 J J n V(k + qy v 2 "J n k 

Pi ~ m t pi - m t 



(6) 

where Y^ va = 2q^g„ a - 2q (2k + q) [10], Cf is the Casimir eigenvalue of the 

fermion representation, pi = Pi + q } and 

Up) = 9>c, I y^JTJ^/ = E """W « 7 » 

is the one loop g.i. quark propagator, derived in [11] 

In principle, this procedure can be generalized to an arbitrary n-point function. In 
particular, the g.i. three and four point functions T^ ua and T^ap have been derived 
in [9] and [12] The Green's functions obtained via the PT, in addition to being g.i., 
are endowed with several characteristic properties. Most noticeably, the gluon n-point 
functions computed thus far [n^ = t^ft, T^q,, Tnua/3 ( n = 2,3,4)] satisfy the following 
simple QED-like Ward identities: 

Qi'r tlva (q 1 ,q2,q3) = t va (q 2 )d~ 1 (q 2 ) - t va (q 3 )d~ 1 (q 3 ) 

= f*»P T l%(<H + ?2, ?3, ?4) + C.p. , 

where d~ 1 (q) = q 2 — fi-(q), f ahc are the structure constants of the gauge group, and the 
abbreviation c.p. in the rhs stands for "cyclic permutations". In addition, the gluon-quark 
vertices [t a ^ of Eq. (6) are ultra-violet finite. 

Regardless of any such properties, however, an ambiguity is associated with the con- 
struction of Green's functions via the PT. It is obvious for instance that, after a g.i. gluon 
self-energy and gluon-quark vertex has been constructed via the PT, one still has the 
freedom to add an arbitrary term of the form (q 2 g iil , — qnqv)f(q 2 ) to the self-energy, and 
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subtract it from the vertex. As long as the function f(q 2 ) is g.L, such an operation satisfies 
the criterion of individual gauge invariance for the self-energy and vertex, respects their 
Ward identities, and preserves the uniqueness of the S-matrix. It is therefore desirable to 
have a physical prescription which eliminates this ambiguity. In this paper we propose an 
alternative, physically motivated prescription for extracting g.i. sub-amplitudes of the S- 
matrix, that are kinematically akin to self-energies, vertices and boxes. In this framework, 
the effective one-loop gluon self-energy is defined to be the limit of the S-matrix as both 
external fermion masses m\ and m 2 are taken to be larger than any other mass scale in 
the process (they are however comparable to each other, e.g. m\ 1712). To one loop we 
show that this limit coincides with the static quark-antiquark potential for very massive 
quarks [13-14], and is physically equivalent to the PT result of Eq. (5). 

Any reasonable effective propagator should only depend on the momentum transfer 
t but not on kinematical details such as masses or total momentum s of the incoming or 
outgoing particles. Similarly, any viable definition of a Gqiqi vertex should only depend 
on t and the quark mass m; [15] and no other kinematical details. This reasoning can 
obviously be generalized to higher n-point functions. Motivated by these observations, 
we propose to define propagator and vertex-like sub-amplitudes by taking appropriate 
kinematical limits of the S-matrix. 

For the simple case of a four-quark on shell amplitude T (Fig.l) we define the following 
three limits 

L (t) = T(s,i,mi = M,m 2 = M) 

Li(t,mi) = T(s,t,mi,m2 = M) (9) 

L 2 (t,m 2 ) = T(s,t,mi = M,m 2 ) 

where the mass M is assumed to be larger than any other mass scale appearing in the 
process, except for any cutoffs introduced in intermediate calculations in order to regularize 
ultra-violet divergences. Note however that, since the external particles are on shell, s = 
(pi +P2) 2 > (itt-i + nT-2) 2 is also of the order of M 2 , in any of these limits [16]. Each of 
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the above quantities is g.L, since it corresponds to a particular limit of the g.i. S-matrix 
element T. They can be systematically computed by expanding the S-matrix in powers of 
(^f-), where Lio is any of the remaining mass scales. The limits considered above correspond 
to well-defined physical situations. To, for example, is the dominant contribution to the 
S-matrix when the momentum transfer t is considerably smaller than the masses of all the 
scattered particles, e.g. t = —q 2 <C mi,m 2 . 

We can define the following linear combinations: 

f (t) = L 

f i {t,m i ) = {L i -L Q ) (z = 1,2) (10) 
f 3 (s,t,mi,m 2 ) = T(s,t,mi,m 2 ) - L - [(T x - T ) + (T 2 - L )] 

We have thus arrived at a decomposition of the S-matrix into individually g.i. and kine- 
matically distinct sub-amplitudes, which we can identify as effective self-energy To(t), 
vertices Ti(t,mi) and T 2 (t,m 2 ), and boxes Ts(s } t, mi , m 2 ). Clearly, the sum of these 
sub-amplitudes is the original S-matrix, e.g. 

T(s,t,mi,m 2 ) = T (t) + Ti(t,mi) + T 2 (t,m 2 ) + T 3 (s,t,mi,m 2 ) , (11) 

The above decomposition of the S-matrix into individually g.i. and kinematically distinct 
sub-amplitudes, relies on a procedure different from the PT. The question that naturally 
arises is how the individual terms of Eq. (2) and Eq. (11) are related. As we will show by 
an explicit one-loop calculation, To of Eq. (2) and To of Eq. (10) are related as follows: 

fo(t) = fo(t) + /u l7 ^ 1 (4)[CV](4)«27^2 • (12) 

q^ 

In Eq. (12) C is a g.i. finite numerical constant. Thus the g.i. self-energy H^ u (q) extracted 
from To(t) , and the To obtained form Eq. (2) satisfy: 

n F ( ? ) = n F ( g ) + Cf F (13) 
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Clearly, the term proportional to C in the r.h.s. of Eq. (13) can be removed by a finite 
counterterm, or, equivalently, absorbed in the final normalization of the S-matrix. 

To compute the leading one-loop contribution to To(t) (or equivalently Lq) we evaluate 
the S-matrix in the limit m\ m,2 — > M, where M ^> —q 2 . [17] For simplicity we consider 
elastic scattering, so that q 2 < 0. As a consequence, there are no imaginary parts in the 
Feynman graphs. We define the Euclidean momentum Q 2 = —q 2 > 0. Throughout the 
calculation we use dimensional regularization, where the UV cutoff is set by the usual pole 
^. In addition, the t'Hooft mass \i has to be introduced. The infrared divergences are 
regulated by introducing an infrared gluon mass A in the intermediate calculations. [18] 

We then compute all one-loop Feynman graphs contributing to the process, neglecting 
terms proportional to any of the ratios (^), an d (^), (or higher powers of such 

ratios), and retaining only logarithmic and constant terms. We emphasize that the above 
expansion is carried out after the integration over the loop momenta has been performed 
in dimensional regularization. Effectively this means that M is always much smaller than 
the cutoff A [e.g. ^ — > hi(^) ^> ln(^j-)]. In this calculation all choices for £ are equivalent, 
since the S-matrix element is ^-independent; we choose £ = 1 for convenience. 

The most involved part of the calculation are the box diagrams. It is important to 
recognize that both the direct and the crossed graph must be appropriately combined 
in order to obtain the correct color structure. It is also interesting to notice that the 
expressions that survive the large M limit are of non-Abelian nature only, namely pro- 
portional to c a . If we call Bdir the total contribution of the direct graph and B cr the 
respective contribution from the crossed, we have that Bdir = {RaRb)i{RaRb)2^dir an d 
B cr = (R a R} ) ) 1 (Ri ) R a ) 2 S cr where Sdir and S cr are the remainders of the boxes, after the 
color structure has been factored out. The important step is to show that in the large M 
limit we have Sdir = —S cr - Thus, the total box contribution B becomes: 

B = (R a Rb) 1 [R a -, Rb] 2 Sdir 

_ 1 (14) 

— r ) C a(Rc)l(Rc)2^dir 
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The result for the individual Feynman graphs are: (we omit external spinors and an overall 

2 

factor of 



[(«)] = 4i =1) 



[(b) + (6) 
[(c) + (c) 
[(d) + (d) 
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mirror \ 



16vr 2 
-.2 



3CW + 4 + 31n(-^) 



9,1V + ••• 



] =77^/-^) 



16vr 2 



A 

C^ + 4 + 21n(-) + 31n(-^) 



16vr 2 



A 

CW + 4 + 21n(-) + 31n(-^) 



(15) 



g 2 ^ 2 

[(e) + (e) crossed )] =2y^-jc a ln( — )g tlu + ... 

where the ellipsis denote terms of order 0( jj) or higher. We notice that the sum of all 

2 o 2 

vertex and box graphs listed in Eq. (15) is equal to 2 \^ [Cjjv ~ l 1 ^^)]? which is, up to a 
physically irrelevant constant, the pinch contribution to the self-energy, given in Eq. (5). 
The total contribution to To reads: 



T = g 2 u 1 -/ tl u 1 (^-) 



^ + 2(0)^ " ln (-#)] V(i) 



g 2 u 1 j^u 1 (^-) 



nr i} + n p (?)-4(f^)v 



(-7)^27^2 



(16) 



1 g 2 c a 1 

T - g 2 uil^ui(-^)[^(j^)^v](-^ )u2l v u 2 



z 

which is the advertised result in the first line of Eq. (12), with C = — 4( \^ )■ The first 
relation of Eq. (13) follows immediately from Eq. (16), namely 



Adding the tree-level contribution to To of Eq. (16), and using the standard result 



(17) 



n (?=l) = 9 C a 

»" 16vr 2 



- 3 [Cuv-H^)] + - 



h {IV 



(18) 
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together with Eq. (5), we find that To is identical to the Fourier transform of the un- 
renormalized one-loop static potential V(Q 2 ) for a heavy quark-antiquark system ( [13]), 
namely (we omit the external spinors), 



where the factor (^f ) 1 g^ r 2 is the coefficient bo in front of —g 3 in the one loop fi function. The 
contribution To of Eq. (16) or Eq. (19) to the S-matrix is infrared finite. The decomposition 
of Eq. (9), together with Eq. (16), implies that the PT result for the self-energy gives the 
dominant contribution to the physical S-matrix, when the scattered particles are heavy 
compared to other mass scales. We have thus arrived at a physical interpretation of this PT 
sub-amplitude. In that sense, the mathematical ambiguity in defining a g.i. propagator- 
like sub-amplitude of the S-matrix, which we discussed previously, can be eliminated by 
imposing a physically motivated boundary condition, i.e. that the effective self-energy 
should reproduce the S-matrix for the scattering of sufficiently heavy external quarks. 
Since perturbation theory in QCD is reliable only for momentum transfers beyond a few 
GeV, in practice this sub-amplitude will provide the dominant contribution to the S-matrix 
only for top and bottom scattering. Nevertheless, it makes sense to define this g.i. sub- 
amplitude also for considerably lighter systems, although in such a case it will generally 
not give the dominant contribution to the S-matrix. 

It would clearly be of interest to extend this analysis to the vertex-like sub-amplitudes. 
Of course, in a theory with massless gauge bosons such sub-amplitudes are in general 
infrared divergent; they can therefore not be directly related to a physical process, without 
including bremsstrahlung. One could nevertheless compare the g.i. vertex-like amplitudes 
Ti and Ti, i = 1,2, of the two schemes, as long as the infra-red singularities are regulated 
in a gauge invariant manner, such as dimensional infra-red regularization [19-20]. This 
goes however beyond the scope of the present communication. 



To = V(Q 2 ) 




a 



1 + 



9 2 Cg 

16vr 2 




(19) 
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The previous considerations can be generalized to the case of multi-quark scattering. 
In particular, from a 2n-quark amplitude one can define a g.i. gluon n-point functions 
r^(<Zi , q n ) with all incoming momenta i = l,..n off-shell. To that end one has to 
consider the limit of the amplitude as all external fermion masses become large (m; — > M, 
i = l,..,n). It would be very interesting to determine if the g.i. n-point functions so 
obtained are physically equivalent to those obtained with the PT, especially for n = 3,4. 
It would be also interesting to generalize the previous arguments to the case of theories with 
spontaneous symmetry breaking in general, and the electro-weak sector of the Standard 
Model, in particular. 
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3. Figure Caption 

Graphs (a)-(e) are some of the QCD contributions to the S-matrix for four-fermion 
processes. Graphs (f), (g), (h), and (i) are the pinch parts of (b), (c), (d), and (e) respec- 
tively. When added to the usual self-energy graphs (a), they give rise to the ^-independent 
amplitude To(t). The mirror image graphs corresponding to (b), (c) and (d), as well as the 
crossed box graph are not shown. Graph (a) contains contributions from virtual fermions, 
gluons and Faddeev-Popov ghosts. Solid (wavy) lines represent fermions (gluons). 
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